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Electromagnetic fed intensity calculation apparatus and method 



(57) In calculating an electromagnetic field intensity 
by a moment method, real and imaginary parts of a mu- 
tual impedance, mutual admittance, and mutual reac- 
tion between elements and of the elements themselves 
are approximated by a sum of terms each with a fre- 
quency raised to a power. The mutual impedances, etc. 
at a plurality of sample frequencies are calculated using 
known mathematical equations, and simultaneous 
equations obtained by substituting the sample frequen- 
cies and the mutual impedances at the respective sam- 
ple frequencies into an approximate equation are solved 
to calculate approximation coefficients for the approxi-. 
' mate equation. Using the thus obtained approximate 
equations, the mutual impedance, etc. at a designated- 
frequency are calculated, and substituted into simulta- 
neous equations describing a boundary condition for 
. each element; by solving the simultaneous equations, 
the electric current, equivalent electric current, and 
equivalent magnetic current of each element are calcu- 
lated, from which the electromagnetic field intensity is 
calculated. 
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Description 

The present invention relates to an electromagnetic field intensity calculation apparatus and an electromagnetic 
field intensity calculation method. More particularly, the invention relates to an electromagnetic field intensity calculation 
s apparatus and an electromagnetic field intensity calculation method that calculate the electromagnetic field intensity 
radiated from an electric circuit device at high speed while maintaining high accuracy. 

Electromagnetic waves radiated from electric circuit devices interfere with other electromagnetic waves such as 
those of television and radio broadcasting, and strict regulations have come to be imposed on such radiation in various 
countries of the world. Examples of standards defining such regulations include the VCCI standards of Japan, the FCC 
w standards of the U.S.A., and the VDE standards of Germany. 

Various techniques such as shielding techniques and filtering techniques are used to meet these electromagnetic 
wave standards. When employing such techniques, it is necessary to carry out a simulation to quantitatively determine 
how much of the energy of the radiation the technique concerned can shield. A need therefore arises for an electro- 
magnetic field intensity calculation apparatus that can simulate electromagnetic field intensities radiated from electric 
J5 circuit devices with high accuracy. 

The strength of an electromagnetic field around an object of an arbitrary shape can be easily calculated using 
known theoretical equations, if the electric current flowing in each part of the object is known. Theoretically, this current 
value can be obtained by solving Maxwell's electromagnetic wave equations under given boundary conditions. How- 
ever, analytical solutions for an object of an arbitrary shape under complex boundary conditions are not available yet. ' 
20 The methods, for obtaining the electric current, which are used in electromagnetic field intensity calculation appa- 

ratus today, all provide approximate solutions. At present, three kinds of approximate solution methods are known: the 
infinitesimal loop antenna approximation method, the distributed constant transmission line approximation method, 
and the moment method. 

The infinitesimal loop antenna approximation method is a method in which the conductor connecting between a 

25 wave source circuit and a load circuit is treated as a loop antenna and, by assuming the electric current on the loop to 
be flat, the electric current is determined using a calculation technique of a lumped constant circuit. Calculation by the 
infinitesimal loop antenna approximation method is the simplest. However, since the calculation accuracy drops sig- 
nificantly under conditions in which the loop dimension becomes significant compared to the wavelength of electro- 
magnetic radiation, this method is of very little use in practice. 

30 The distributed constant transmission line approximation method is a method in which the electric current is ob- 

tained by applying a distributed constant transmission line equation to an object that can be approximated as a one- 
dimensional structure. Calculation by the distributed constant transmission line approximation method is relatively 
simple, with the required computation time and memory capacity only increasing almost proportionally to the number 
of elements to be analyzed; furthermore, phenomena such as reflection and resonance of the transmission line can 

35 also be analyzed. An analysis can therefore be done with high accuracy and at high speed for objects for which one 
dimensional approximations are possible. However, the distributed constant transmission line approximation method 
has a problem in that an analysis cannot be done for objects that cannot be approximated to one-dimensional structures. 

The moment method is a method of solving an integral equation derived from Maxwell's electromagnetic wave 
equations, and can handle three-dimensional objects of any shape. More specifically, the moment method divides an 

40 object into small elements (wires or surface patches or the like) and calculates the electric current flowing on each 
segment, thereby computing the electromagnetic field intensity. Since the moment method can handle three-dimen- 
sional objects of any shape, the design of electromagnetic field intensity calculation apparatus predominantly employs 
a configuration in which the strength of an electromagnetic field radiated from an electrical circuit device is calculated 
using the moment method. 

45 Reference literature on the moment method includes the following. 

H. N. Wang, J. H. Richmond and M. C. Gilreath: "Sinusoidal Reaction Formulation for Radiation and Scattering 
from Conducting Surface, ' IEEE TRANSACTIONS ANTENNAS PROPAGATION, AP-23, pp. 376 - 382, 1975. 

In the moment method, as described for example in Japanese Patent Unexamined Publication Nos. 7-234890 and 
7-302278 (U.S. Ser. No. 432,261) by the present inventor et al., a conductor constituting an electric circuit device, 
. so including a housing, cable, etc., is divided into small elements such as wires or surface patches. The mutual impedance 
between the elements and the self-impedance of each individual element (collectively called the mutual impedance 
Zy) are calculated using known mathematical equations from frequency values and the geometric data of the elements 
concerned, and the value of the mutual impedance Zjj thus obtained is substituted into simultaneous equations de- 
scribing a boundary condition for each element. By solving the simultaneous equations, the current flowing on each 

55 element is determined, from which the electromagnetic field intensity is calculated. When it is desired to increase the 
accuracy by taking into account the scattering of the electromagnetic field caused by a dielectric contained in the 
electric circuit device, the dielectric also is divided into small elements, and the mutual admittance and mutual 
reaction Bj, are calculated for each element in addition to the mutual impedance Z^; then, simultaneous equations 



2 



EP 0 778 533 A2 



containing a boundary condition 1or each element of the dielectric are solved. 

In specifications defining EMC electromagnetic wave regulations, allowable values are specified for the regulated 
frequency range. For example, the VCCI standards, defining EMC electromagnetic wave regulations in Japan, stipulate 
the regulated frequency range of 30 MHz to 1 GHz and specify allowable values for that frequency range. Consider a 
s transmission line on a printed circuit board as an example of an electric circuit device subject to the electromagnets 

- wave regulations. The transmission line sends an output of a driver (output circuit) to a receiver (receiving circuit). The 
output of the driver has a pulse-like voltage waveform. The components actually radiated from the transmission line, 
therefore, contain the fundamental frequency (f 0 ) component of the clock, plus its harmonic components (f 1t 

f ) whose frequencies are integral multiples of the fundamental frequency. As a result, in an electromagnetic field 

10 intensity calculation apparatus, calculations for the simulation of electromagnetic field intensities must be performed 
on all the harmonic components generated from the transmission line and falling within the regulated frequency range. 
An electric field spectrum is obtained as a result of such simulation. The resulting electric field spectrum is a set of 
electric field intensities calculated at the respective frequencies by the electromagnetic field intensity calculation ap- 
paratus using the moment method. Accordingly, to obtain the electric field spectrum, the electromagnetic field intensity 

is calculation by the moment method has to be performed for each frequency (each of the fundamental frequency and/ 
or harmonic frequencies). That is, the calculation process has to be repeated for each frequency. This has lead to the 
problem that it takes a very long time for the prior art electromagnetic field intensity calculation apparatus to determine 
through simulation whether the EMC electromagnetic wave specifications are satisfied or not. For the calculation of 
mutual impedances, in particular, there has been the problem that the calculation takes an extremely long time because 

20 the amount of calculation is enormous. 

Thus the prior art electromagnetic field intensity calculation apparatus employs a method that calculates the mutual 
■ impedance for 'each" frequency by carrying out known processing commonly used in the moment method for the 
calculation o1 the mutual impedance. The prbr art electromagnetic field intensity calculation apparatus therefore has 
had the problem that processing involving a considerable amount of calculation has to be performed for each frequency 

25 This problem becomes more serious when the mutual admittance or mutual reaction has to be obtained in addition to 

the mutual impedance. . , '*u 

An example of calculation time will be given below. The calculation time depends much on the shape of a three- 
dimensional object to be simulated. To give a rough value, when a three-dimensional object to be simulated was divided 
into elements consisting of 1 000 surface patches, it took several hours to calculate the electromagnetic field intensity 
30 for each frequency. More specifically, it took several hours to calculate the mutual impedance, several minutes to 

- calculate the simultaneous equations, and several minutes to calculate the electric field or magnetic field. 

It is accordingly desirable to provide an apparatus for calculating an electromagnetic field intensity and a method 
of calculating the same, capable of calculating the strength of an electromagnetic field radiated from an electric circuit 
device at high speed while maintaining high accuracy. ' 

35 According to an embodiment of the present invention, there is provided a method of calculating an electromagnetic 

. field intensity, comprising the steps of: a) calculating approximation coefficients for each of approximate equations that 
respectively yield approximate values of electromagnetic characteristic values, between a plurality of elements con- 
stituting an electric circuit device and of the elements themselves, at an arbitrary frequency; and b) calculating the 
app roximate values of the elect romagnetic characteristic values between the elements and of the elements themselves 

40 • at a designated frequency in accordance with the approximate equations having the approximation coefficients calcu- 
lated in step a), thereby making it possible to evaluate the electromagnetic field intensity for the designated frequency. 

According to an embodiment of the present invention, there is also provided an apparatus for calculating an elec- 
tromagnetic field intensity, comprising: means for calculating approximation coefficients for each of approximate equa- 
tions that respectively yield approximate values of electromagnetic characteristic values, between a plurality of ele- 

45 ments constituting an electric circuit device and of the elements themselves, at an arbitrary frequency; and means for 
calculating the approximate values of the electromagnetic characteristic values between the elements and of the ele- 
ments themselves at a designated frequency in accordance with the approximate equations having the approximation 
coefficients calculated by the approximation coefficients calculating means, thereby making it possible to evaluate the 
electromagnetic field intensity for the designated frequency. 

so According to an embodiment of the present invention, there is also provided a program storage device readable 

by a machine, tangibly embodying a program of instructions executable by the machine to perform method steps for 
calculating an electromagnetic field intensity, the method steps comprising: a) calculating approximation coefficients 
for each of approximate equations that respectively yield approximate values of electromagnetic characteristic values, 
between a plurality of elements constituting an electric circuit device and of the elements themselves, at an arbitrary 

55 frequency and b) calculating the approximate values of the electromagnetic characteristic values between the elements 
and of the elements themselves at a designated frequency in accordance with the approximate equations having the 
approximation coefficients calculated in step a), thereby making it possible to evaluate the electromagnetic field inten- 
sity for the designated frequency. 
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Reference will now be made, by way of example, to the accompanying drawings, in which:- ~ 

Figure 1 is a block diagram of an electromagnetic field intensity calculation apparatus according to the present 
invention; 

s Figure 2 is a diagram showing an example of a frequency spectrum of electromagnetic field intensities, output 

from the electromagnetic field intensity calculation apparatus; 

Figure 3 is a flowchart illustrating the operation of the electromagnetic field intensity calculation apparatus; 
Figure 4 is a flowchart for the calculation of approximation coefficients by LU decomposition; 
Figure 5 is a diagram showing simultaneous equations accounting for the effect of scattering caused by dielectrics; 
10 Figure 6 is a diagram showing simuttaneous equations further accounting for the presence of a core inserted in a 

cable; 

Figure 7 is a diagram showing a matrix equation for a physical object; 
Figure 8 is a diagram showing a matrix equation for a mirror image; 
Figure 9 is a diagram showing simultaneous equations accounting for reflected waves; 
is Figure 1 0 is a diagram showing sample frequencies for approximate calculations of impedance and a frequency 

for approximate calculation; 

Figure 1 1 is a diagram showing monopoles assumed in a comparison test; 

Figures 12 to 18 are diagram showing results of the comparison test; and 

Figure 1 9 is a diagram showing monopoles used in the proof of an approximate equation. 

20 

Figure 1 shows the configuration of an electromagnetic field intensity calculation apparatus 20. The electromagnetic 
field intensity calculation apparatus 20 obtains mutual impedances, etc. using approximate equations, and by using 
the obtained results, solves simultaneous equations of the moment method and thereby computes the strength of an 
electromagnetic field radiated from an electric circuit device to be analyzed. The electromagnetic field intensity calcu- 

25 lation apparatus 20 comprises a CPU, memory, and electromagnetic field intensity calculation programs stored in the 
memory. The electromagnetic field intensity calculation programs may be provided on a known storage medium such 
as a flexible disk or a CD-ROM. 

When the electromagnetic field intensity calculation apparatus 20 is activated, an input processing section 1 reads 
structural information stored in an input file located outside the figure. The structural information is needed to apply 

30 the moment method to the electric circuit device for which a simulation is to be performed, and is used to divide the 
electric circuit device into a plurality of elements (surface patches or wires). 

Various data necessary for electromagnetic field intensity calculations are input to the input processing section 1 
from outside the electromagnetic field intensity calculation apparatus 20. When a frequency range and frequencies to 
be simulated (analyzed) are input, these are stored in a memory located outside the figure by the input processing 

35 section T. 

The frequency range defines a simulation range. For example, when carrying out an analysis to determine whether 
the electric circuit device to be analyzed satisfies the VCCI specifications, the frequency range is set equal to the 
control frequency range of 30 MHz to 1 GHz. The frequency range can also be set to define any desired simulation 
range, regardless of the EMC electromagnetic wave regulation specifications. 

40 The frequencies to be simulated are the frequencies at which the electric field and magnetic field intensities are 

to be calculated. For example, frequencies falling within the input frequency range and equal to integral multiples of 
the fundamental clock frequency of the electric circuit device to be analyzed are input as the frequencies to be simulated. 

Usually, there are many frequencies to be simulated within the input frequency range. First, the fundamental clock 
frequency of the electric circuit device to be analyzed, for example, is input to the input processing section 1. Then, 

<5 based on this fundamental clock frequency, the input processing section 1 calculates all the frequencies that are equal 
to integral multiples of the fundamental clock frequency and that fall within the frequency range. These frequencies 
are the frequencies to be simulated, and the mutual impedance at each frequency is obtained using an approximate 
equation of the present invention. The input processing section 1 stores the thus calculated frequencies to be simulated 
into the memory located outside the figure. 

50 Alternatively, the frequencies to be simulated may be designated specifically and input directly to the input processing 
section 1, as previously described. 

A sample designation processing section 3, a sample calculation processing section 4, and an approximation 
coefficient calculation processing section 5 together constitute an approximation coefficient calculation section 2. The 
approximation coefficient calculation section 2 contains a data file 6 which consists of an impedance file 7 for storing 

55 mutual impedances and a coefficient file 8 for storing approximation coefficients. 

The sample designation processing section 3 determines sample frequencies based on the frequency range and 
the number of samples, n. For this purpose, the sample designation processing section 3 refers to the frequency range 
stored in the memory. The sample frequencies are samples used for the calculation of approximation coefficients. A 
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oluralitv of such sample frequencies (n sample frequencies) are designated. The value n is determined according to 
Curacy of the approximate equation used by an approximate calculation processing secfoo 9 H will be described 
ater More n is set to a value not smaller than the number of terms in the real and the imagery part of 
■ he approximate equation, and preferably to a value equal to the number of terms. The following description deals w, h 
s TLTe in" n ? s chosen to be equal to the number of terms in the real and the imaginary part o the approximate 
equate The sample designation pressing sectbn 3 stores the thus determined sample frequences ,n the memory 

'° Ca Tf e^mJe^reS^ies are designated a, substantially equispaced frequency interval with spacing .chosen to 
• be equal to the width of the frequency range divided by (n-1 ). For example, when the frequency range is from 30 UHz 
« to 1 GHz and n = 5. frve sample frequencies, t,, . f^ fc. U and U are chosen as specie numerical values. In th,s 
example, the sample frequencies. f s1 , i S2 . f*. U and U are chosen as 30 MHz, 250 MHz. 500 MHz and 750 MHz 
• resSely. The lowest sample frequency , sl is not 0 MHz. which would be the value if chosen at 

intervals, but is chosen as 30 MHz, thesmallestvalue of thecontrol frequency range. Howa^equispacmg the sample 
frequenc es is not strictly demanded. After all. it is only necessary to construct simultaneous equations tor obtaining 
« approximation coefficients, and if the spacing is somewhat unequal, it will not present a problem ,n P«^B- 

The sample calculation processing section 4 calculates electromagnetic characteristic values, such as the mutual 
impedance mutual admittance, and mutual reaction, at n sample frequencies. For this purpose, the sample ca cu a on 
pTo P cessings™4reters,othensar^^ 
' 'he electromagnetic field intensity without considering the effect of the electromagnetic field sc «^ 
20 dielectrics, mlal impedances on* for the elements of the conductor (metal are calculated and when the effec ; s 
to be considered, muiual impedances, mutual admittances, and mutual reactions are calculated for all the elements 
needing the elements of the dielectrics. The mutual impedance, etc. are calcuteted from the 9-™rtn^tatt £ the 
elements concerned and frequency values, without using approximate equations, but by using known mathematical 
"ns commonly used in the moment method. The values of the calculated mutual impedance, etc. are therefore 

" When only the mutual impedance Z, is to be calculated, a number, n. of mutual impedances ^ •*^* a ™"*"' 

■ of sample frequencies, are calculated tor each pair of elements determined by the values of , ^J™*™*^ 
culation processing section 4 stores the calculated mutual impedance values .n the impedance file 7. When the mutual 
admittance and mutual reaction are calculated, these also are stored in the file 7. = , imr , oHa n,.« 

30 When n = 5 with five sample frequencies designated as f,, , f s2 . W W and W for example, the mu ua -mpsdances 

Z- , to Zr rf are obtained for the respective sample frequencies of f s1 to W Thus, specific values of mutual impedances 
Z.'' to Z^ are obtained for each pair of elements determined by the values of i and j. These mutual .mpedances 
Zj ii5l to Zjj iS5 are expressed as follows: 
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Z ij^3 = a q*3 + iPii53 



Z ii.s4 = a ii.s4 +iP ii*' 



Z ij.sS = a ii.sS + iPii^S 



so . where a, is the value of the real part of the mutual impedance Z* ft, is the value of the imaginary part of the mutual 
££2£z, and I and j, the subscripts on a and p. each take a value between t and m ™«^toenunte« 
■ modes of the moment method. Further, j prefixed to p is the imaginary unit and .s distinguished from the , w itten 
«Uo p. Thus, five sets of real part a, and imaginary part Pij are respectivefy obtained as specific numeric ^ 
Usinq the n sets of mutual impedances, etc. calcubted by the sample calculation processing section 4 and the n 
ss sample frequencies designated by the sample designation processing section 3, the approximation coefficient calcu- 
aZ processhg section 5 calculates approximation coefficients in accordance with the ^^^f^^ 
present invention. For this purpose, the approximation coefficient calculation processing section 5 refers to the n sets 
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of mutual impedances, etc. and the n sample frequencies stored in the memory. 

According to the present invention, the approximate equation for the mutual impedance Zg is 

5 Z- = exp(-jkr)[a 0 + a, f 2 + a 2 f 4 + a 3 f 6 + a 4 f 8 +... + j^bj/f + b 1 f + b 2 f 3 + b 3 f 5 + b 4 f 7 + ...)] 

where 1 is the frequency, Bq, a,, ... are approximation coefficients for the real part, and b 0 , b 1 , ... are approximation 
coefficients for the imaginary part. Further, j is the imaginary unit, k is a value determined by k = 2rcf/c, and r is the 
distance between the centers of two elements. Needless to say, the term "a^ is actually the term *a 0 xi° = a^, and 
10 the term "b(/f" is actually the term "b 0 x f 1 ". Accordingly, the real part is a sum of the terms with the frequency raised 
to an even-numbered power (0, 2, 4, 6, 8, ...) and multiplied by the respective approximation coefficients, and the 
imaginary part is a sum of the terms with the frequency raised to an odd-numbered power (-1,1, 3, 5, 7, ...) and 
multiplied by the respective approximation coefficients. 

The approximate equation for the mutual admittance is 

15 

Y- = exp(-jkr)[a 0 + a^ 2 + a 2 f 4 + a 3 f 6 + a 4 f 8 + ... + \{bji + b n f + b 2 f 3 + b 3 f 5 + b 4 f 7 + ...)] 

which is of the same form as that for the mutual impedance. 
20 The approximate equation for the mutual reaction Bjj is 

Bjj = exp(-jkr}[c 0 + c.f 2 + c 2 f 4 + c 3 f 6 + c 4 f 8 + ... +. j(d 0 f + d/ + d 2 f 5 + d 3 f 7 + d 4 f 9 + ...)] 

25 where c 0 , c 1( ... are approximation coefficients for the real part, and d 0 , d v ... are approximation coefficients for the 
imaginary part. This equation is similar in form to that for the mutual impedance, the difference being that the first term 
of the imaginary part begins with the 1 st power of the frequency f instead of the -1 st power of the frequency f ; 

How the mutual impedance, etc. can be approximated using the above approximate equations will be explained 
later. 

so The approximate equations are given in advance to the approximation coefficient calculation processing section 

5 in the form of a program. At the same time, the number, n, of terms in the real and the imaginary part {hence the 
number, n, of samples) is determined. For example, when n = 5, the approximate equation for the mutual impedance 
Zjj is given as 

35 Zjj = exp(-jkr)[a 0 + a, f 2 + a/ + a 3 f 6 + a 4 f 8 + j(b 0 /f + b 1 f + b 2 f 3 + b 3 f 5 + b/)] 

Likewise, the approximate equation for the mutual admittance Y rj is given as 
40 Y- = exp(-jkr}[a 0 + a 1 f 2 + a 2 f 4 + a 3 f 5 + a 4 f 8 + j(b 0 /f + b 1 f + b 2 f 3 + b 3 f 5 + b 4 f 7 )] 

The approximate equation for the mutual reaction Bjj is given as 
45 ' B- = exp(-jkr)[c 0 + c, f 2 + c/ + c 3 f 6 + c 4 f 8 + j(d 0 f + d/ + d 2 f 5 + d/ + d 4 f 9 )] 

Provisions may be made so that the number, n, of terms in the approximate equations (hence the number, n, of 
samples) can be specified or changed by an input from the input processing section 1. It is generally believed that the 

so accuracy of approximation depends on the number of terms, n, in an approximate equation. Therefore, the form of the 
approximate equation, that is, the number of terms, n : in the approximate equation, may be changed according to the 
situation; that is, when the accuracy of approximation is not an important concern, the number of terms, n, can be 
reduced to speed up the calculation of electric fields, etc., and when the processing time is not an important concern, 
the number of terms, n, is increased so that electric fields, etc. can be obtained with higher accuracy. 

ss For the calculation of the approximation coefficients, the approximation coefficient calculation processing section 

5 reads out the n sets of mutual impedances, etc. stored in the impedance file 7, but after they have been read out, 
the impedance file 7 is no longer necessary. Accordingly, the approximation coefficient calculation processing section 
5 stores the calculated approximation coefficients in the file used as the impedance file 7. More specifically, the im- 
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pedance file 7 is now used as the coefficient file B. That is, the same file is used as the impedance file 7 at one time 
' and as the coefficient file 8 at other times. Since the data amount of the mutual impedance, etc. and approximation 
coefficients to be stored in a file is large, memory resources can be used effectively by sharing the same file in this 

s man For example, when n=5 has been given and five sets of mutual impedances to Z jjs5 have been obtained for 
five sample frequencies f s1 to f^, the approximation coefficient calculation processing section 5 constructs five simul- 
taneous equations for the real part by substituting the five values of the real part and the five sample frequencies 
into the real part of the approximate equation. These simultaneous equations are for (the real part a\ of) the mutual 
impedance Zjj between the two elements determined by the values of the subscripts i and j. The simultaneous equations 

io are expressed as 

ct'^ =a 0 + a 1 f s1 2 + a 2 f s1 4 +a 3 f sl 6 + a 4 f 3l 8 (2) 



15 



20 



a 'ij.s2 = 3 0 + 3 1 f s2 2 + a 2 f s2 4 + a 3 f s2 6 + a 4 f s2 

a 'ij,s3 = a 0 + a i f s3 2 + 3 2 f s3 4 + 3 3 f s3 6 + ' a 4** 
2 , 4 . 6 „ , 8 



In the above simultaneous equations, the values of the real parts a' ?iSl to a' ijs5 of the mutual impedances are calculated 
from the following equations by using the Zy iSl to Z ij s5 calcu lated in the sample calculation processing section 4. . 

30 - . . . Z^'expH^^ (3) 

Z iis2 /exp[-j(27cf s2 /c)r]= a'jj^+iP'jj^ 

35 

Z-j s4 /exp[-i(27cf s4 /c)r] = a' ijs4 + tf' ijfS4 

40 " 

z r h s5 /ex PH( 2 ^W c ) r ] = a, ij,s5 + iP'ii^s 

45 Here, the values of the sample frequencies f 8l to f^ are already designated. Further, the distance r between the centers 
of the two elements determined by the values of i and j is a known value. Therefore, the value of the left-hand side of 
each equation (3) can be calculated, and by substituting the resulting values of the real parts a' jjs1 to a' ji sS in equations 
(2), the five simultaneous equations are solved for the five unknowns slq to a 4 , to obtain the values of the approximation 
coefficients a 0 to a 4 . These approximation coefficients a 0 to.a 4 are used in the approximate equation expressing (the 

. so real part of) the mutual impedance Zy, at designated frequency t, between the two elements determined by the values 
of i and j. 

In like manner, the approximation coefficient calculation processing section 5 constructs five simultaneous equa- 
tions for the imaginary part by substituting the five values of the imaginary part P' g and the five sample frequencies into 
the imaginary part of the approximate equation. These simultaneous equations are for the imaginary part $\ of the 
55 mutual impedance Z {] between the two elements determined by the values of i and j. The simultaneous equations are 
expressed as 
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Msi + Vsi + b 2 f s1 3 + b 3 f s1 5 + b 4 f sl 7 (4) 

+ b 1 f s2 + b 2 f s2 3 + b 3 f s2 5 + b 4*wJ 

b 0^ + b 1 f s3 + b 2 f s3 3 + b 3 f s3 S + b JsJ 

b 0 /J s4 + b 1 f s4 + + b 3 f s4 5 + b 4*s4 7 

MsS + b 1 f s5 + b 2 f s5 3 + ^sS 5 + b 4 f sS ? 
15 

These simultaneous equations are solved to obtain the values of the approximation coefficients b 0 to b 4 . These ap- 
proximation coefficients b 0 to b 4 are used in the approximate equation expressing the imaginary part of the mutual 
impedance 2^, at designated frequency f, between the two elements determined by the values of i and j. The simulta- 
neous equations (2) and (4) can be solved, for example, by the known method of Gaussian elimination, but the solutions 
20 can be obtained much faster by using the method of LU decomposition. Details of the LU decomposition method will 
be described later. 

The processing in the approximation coefficient calculation section 2 has been descrtoed above. The processing 
performed in the approximation coefficient calculation section 2 is to calculate the approximation coefficients used in 
the approximate equations for the mutual impedance, etc. and thereby complete the approximate equations. This 

25 processing is actually preprocessing preparatory to the calculation of the mutual impedance, etc. by approximation. 
This processing is performed only once for one set of n sample frequencies. In the approximation coefficient calculation 
section 2, the processing for the calculation of the mutual impedance, etc. is performed by using known methods 
without using approximate equations. However, since this processing is performed only for a limited number of mutual 
impedances, etc. for each attention element pair, this processing does not have a significant influence on the total 

30 processing time. 

On the other hand, the actual calculation process for the mutual impedance, etc. is performed in the approximate 
calculation processing section 9. The calculation process uses approximate equations of simple form and performs 
calculations by approximation. Accordingly, the time required to calculate the mutual impedance at one frequency is 
extremely short. This calculation processing is repeated for each designated frequency within the specified frequency 

35 range (for example, 98 frequencies as will be described later). However, since the processing time for one frequency 
is extremely short, the total processing time can be reduced drastically. 

Using the plurality of approximation coefficients' calculated by the approximation coefficient calculation processing 
section 5, the approximate calculation processing section 9 calculates the approximate values of the mutual impedance, 
etc. for the designated frequencies within the frequency range in accordance with the previously given approximate 

40 equations. For this purpose, the approximate calculation processing section 9 refers to the coefficient file 8. In the 
previously given approximate equations, values other than the frequency f are already determined. Therefore, by sub- 
stituting a designated frequency into the approximate equations, the mutual impedance, etc. at that frequency can be 
calculated. The approximate calculation processing section 9 stores the calculated mutual impedance, etc. at desig- 
nated locations in the memory. 

45 When n = 5, for example, the approximate equation for the mutual impedance 2,-j is given as 

Zjj = exp[-j(2sf/c)r)][a 0 + a n f 2 + a/ + a 3 f 6 + a 4 f 8 + Kb^f + b, f + b 2 f 3 + b 3 f 5 + b 4 f 7 )} 

so |n this approximate equation, the distance r between the centers of two elements and the approximation coefficients 
ao to a 4 and b 0 to b 4 are known values. Then, by substituting frequency f into the approximate equation, the mutual 
impedance at that frequency can be calculated. 

Using this approximate equation, calculation of the mutual impedance is repeated for each designated frequency 
within the frequency range! For example, when the fundamental clock frequency of an electric circuit device is 10 MHz 

55 and the frequency range is from 30 MHz to 1 GHz, the mutual impedance is obtained for each of the harmonic com- 
ponents of the fundamental frequency of 10 MHz. Since a harmonic has a frequency that is an integral multiple of the 
fundamental frequency, the harmonic frequencies falling within the frequency range of 30 MHz to 1 GHz are 30 MHz, 
40 MHz, 50 MHz 990 MHz, and 1 GHz. That is, there are '98 harmonic components within the frequency range. 
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The mutual impedance at each of the 98 harmonic components is calculated using the approximate equation, rather 
than an exact equation as used in the prior art. ._',,« ^ * 

' Using the electromagnetic field intensity calculation apparatus 20, only mutual impedances may be calculated if 
the role of evaluating the electromagnetic field intensity generated by an electric circuit device is left to other means. 
s m other words, the electromagnetic field intensity calculation apparatus 20 can be used as a mutual impedance cal- 
culation apparatus. In that case, the processing in the electromagnetic current calculation processing section 10 and 
electric/magnetic field calculation processing section 11 is omitted, and only the mutual impedances calculated in the 
approximate calculation processing section 9 are output. In like manner, only mutual admittances or mutual reactions 
may be calculated. . . 

io The electromagnetic current calculation processing section 1 0 substitutes the approximate values of the mutual 

impedance calculated by the approximate calculation processing section 9 into simultaneous equations of the moment 
method expressing a boundary condition for each element, and calculates the value of the current flowing on each 
alement by solving the simultaneous equations. For this purpose, the electromagnetic current calculation processing 
section 10 refers to the approximate values of the mutual impedance that the approximate calculation processing 
is section 9 has stored in the memory, The electromagnetic current calculation processing section 1 0 stores the calculated 
current values at designated locations in the memory. 

The simultaneous equations of the moment method for the mutual impedance 2^, when the effect of the dielectrics 
is not taken into account, are given by [2,^ = [Vj], where [ ] represents a vector or matrix. The earlier calculated mutual 
impedances 2j, are substituted into this matrix equation. Since wave source [Vj] is known, the value of the current [IJ 
20 flowing on each element is obtained. The calculation of the current value [IJ is repeated for each of the designated 
frequencies within the frequency range. , 

The electric/magnetic field calculation processing section 11, using the electric current values calculated by the 
electromagnetic current calculation processing section 1 0, performs prescribed calculations to calculate electric field 
E and magnetic field H. For this purpose, the electric/magnetic field calculation processing section 11 refers to the 
25 values that the electromagnetic current calculation processing section 1 0 has stored in the memory. The electric/mag- 
netic field calculation processing section 11 stores the calculated electric field E and magnetic field H in the output file 
12 In this way. the electromagnetic field intensity radiated from the electric circuit device is calculated. 

When taking into account the effect of the electromagnetic field scattering caused by the dielectrics contained in 
the electric circuit device, the electromagnetic current calculation processing section 10 substitutes the values of the 
so ' mutual impedance, mutual admittance, and mutual reaction, calculated by the approximate calculation processing 
section 9 into simultaneous equations of the moment method containing a boundary condition for each element of the 
dielectric and by solving the simultaneous equations, calculates the equivalent electric current and equivalent magnetic 
current flowing on each element of the dielectric together with the electric current flowing on each element of the 
conductor. From the obtained results, the electric/magnetic field calculation processing section 11 calculates the elec- 
35 tromagnetic field intensity. The simultaneous equations used for the calculation will be described later. 

After that, based on the data of the electric field E and magnetic field H stored in the output file 12, an output 
processing section, which is located outside the figure, creates a frequency spectrum for output in the form of a graph 
which plots the frequenby distribution of the electromagnetic field intensity, as shown in Figure 2. Alternatively, the 
output processing section may be constructed to output the frequency distribution of the electromagnetic field intensity 
40 in other form ! based on the data of the electric field E and magnetic field H stored in the output file 1 2. 

Process Flow for Electromagnetic Field Intensity Calculation _ _ - 

Figure 3 shows a process flow for the electromagnetic field intensity calculation performed in the electromagnetic 
45 field intensity calculation apparatus 20. This process flow shows an example where mutual impedances are calculated 
first and then the electromagnetic field intensity is computed using the mutual impedances. Fundamentally the same 
process applies for cases where mutual admittances and mutual reactions are calculated along with mutual imped- 
ances. . 

When the electromagnetic field intensity calculation apparatus 20 is activated, the input processing section 1 reads 
so the structural information. After that, when a frequency range is specified and entered from the outside, the input 
processing section 1 obtains the frequency range (S1 ). Then the sample designation processing sect.on 3 determines 
a plurality of sample frequencies (S2), based on the frequency range and on the number of terms, n, in the real and 
the imaginary part of the approximate equation to be used in the approximate calculation processing section 9. The 
sample calculation processing section 4 examines whether the calculation of the mutual impedance 2» has been com- 
ss pieted for all of the designated sample frequencies (S3). If it has not been completed yet, the sample calculat.on 
processing section 4 calculates the mutual impedance at the next sample frequency remaining unprocessed, by 
performing known calculations used in the moment method (54). The process then returns to step S3. 

When the calculation has been completed for all the sample frequencies, the approximation coefficient calculation 
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processing section 5 constructs simultaneous equations by substituting each sample frequency and its corresponding 
mutual impedance Zq into the real and imaginary parts of the approximate equation to be used in the approximate 
calculation processing section 9, and calculates the approximation coefficients by solving the simultaneous equations 
(S5). This completes the preprocessing, which is followed by a process in which approximate calculation of the mutual 

5 impedance for each frequency and the calculation of the electric current value using the same are repeated to 
calculate the electric and magnetic fields. 

First, the approximate calculation processing section 9 examines whether the calculation of the mutual impedance 
Zjj has been completed for all the frequencies within the specified frequency range (S6). tf completed, the electromag- 
netic field intensity calculation process is terminated. Thereafter, the electric field spectrum is output If not completed 

10 yet, the approximate calculation processing section 9 calculates the approximate value of the mutual impedance 
for one frequency in accordance with the approximate equation (S7). The electromagnetic current calculation process- 
ing section 10 substitutes the approximate value of the mutual impedance Z^ into the simultaneous equations of the 
moment method, and calculates the electric current value by solving the simultaneous equations (SB). The electric/ 
magnetic field calculation processing section 11 examines whether the calculation of the electric current has been 

is completed for all the observation points (S9). If completed: the process returns to step S6. If not completed yet, the 
electric/magnetic field calculation processing section 11 calculates the electric field E and magnetic field H for the next 
observation point remaining unprocessed, and stores the resutts in the output file 12 (S10). After that, the process 
returns to step S9. 

20 Speeding the Calculation of Approximation Coefficients by LU Decomposition Method 

When directly solving the simultaneous equations (2), (4) by Gaussian elimination or the like, since the values of' 
the real parts a'^ to vary with the values of i and j, the process for calculating the approximation coefficients 
to a 4 by solving the simultaneous equations must be repeated the same number of times as the number of possible 

25 combinations of i and j values. On the other hand, when solving the simultaneous equations by the LU decomposition 
method, though the values of the real parts a'jj s1 to a'jj iS5 vary with the values of i and j, part of the process for solving 
the simultaneous equations can be made common to any combination of i and j values, and furthermore, the remaining 
part of the solving process can be made relatively simple. Accordingly, the approximation coefficients can be calculated 
at high speed by repeating only the remaining part of the solving process in accordance with the variation of the i and 

30 ] values. 

The reason that the approximation coefficients for the approximate equations for the mutual impedance Z,-j, etc. 
can be calculated using the LU decomposition method will be explained below. The following explanation deals with 
the calculation of the approximation coefficients ao to a 4 for the real parts a'jj s1 to a',-^ of the mutual impedance Z,-j. 
In equation (2), the real parts a'^ to a' jjs5 varywith the values of i and j. The sample frequencies f sl to f^ are 
35 fixed values. The five approximation coefficients aQ to a 4 are unknowns. Accordingly, the right-hand side of equation 
(2) takes the form of the unknowns, ao to a 4 , respectively multiplied by the coefficients f sl ° to f^ 8 , etc., and the left- 
hand side of equation (2) can be regarded as variables, a' ri sl to a',^, whose values vary with the values of i and j. 

When equation (2) is seen from such a viewpoint, it is found that the right-hand side, excluding the unknowns ao 
to a 4 , does not change even if the values of i and j change. In other words, even when the values of the real parts a'^ s1 
40 t° a 'q ss of the mutual impedance Z,-j change as a result of a change in the values of i and j, the 5X5 coefficient matrix 
consisting of the coefficients f s1 0 to f^ 8 , etc. on the right-hand side remains unchanged. In equation (2), those values, 
that change according to the values of i and j are the unknowns on the right-hand side and the constant term on the 
left-hand side. 

Noting that the coefficient matrix on the right-hand side of equation (2) does not change with the values of i and j, 
45 the LU decomposition method can be applied to the solution of equation (2) as the simultaneous equations. That is, 
according to the LU decomposition method, the 5 X 5 coefficient matrix consisting of the coefficients f s1 0 to f^ 8 , etc. 
. on the right-hand side is decomposed as a product of a lower triangular matrix L, with the leading diagonal elements 
consisting of 1 , and an upper triangular matrix U, as shown in equation (5) below. 

so 
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Here, for the convenience of expressing equation (6) given hereinafter, Jthe coefficient matrix (the 5X5 coefficient 
matrix consisting of the coefficients f sl ° to f^ 8 , etc.) on the left-hand side of equation (5) is set as follows: 



I 3 l z 

*-*z H ^ 

£ S 3 2 £ *3 A ts3^ 

f,« 2 * s < 4 £i4 6 -f-s 

. r»5 £ s5 r s5 r »3 





C,2 


Cu 


Cu.. 


C13 




1 


C21 


C 2 2 


C 2 3 


C 2 A 


c 25 




1 


C31 


C32 


C33 


C 3 4 


c 35 




. 1 


c 4l 


c« 


C 63 


Ca4 


c« 




1 


C 3l 


C 52 


c 5 , 


C 5 4 • 


c 33 




1 



Then, according to the LU decomposition method, the elements in the lower triangular matrix L and upper triangular 
30 matrix U on the right-hand side of equation (5) can be expressed as shown in equation (6) below. Here, U (n , etc. denote 
the elements in the upper triangular matrix U in equation (5), and L^,, etc. represent the elements in the lower triangular 
matrix L in equation (5), where subscript 1n, for example, denotes the n-th column in the first row, and subscript ml 
indicates the first column in the m-th row. 



35 



U 1n = C m 



(n = 1 to 5) 



(6) 



40 



45 



SO 



L m1 = C ml /U ln <m = 21o5) 
U 2 „ = C 2n -L 21 U,„ (n = 2to5) 
L^MC^-L^U^VU^ (m = 3.4,5) 
U3„ = C3„-L 31 U ln -^U 2n (0 = 3.4.5) 
L m3 = (C m3 • L m , U, 3 - L m2 U 23 yU33 (m = 4, 5) 
U 4 „ = C 4 „ " L 4 , U 1n • L^Ua, - L 43 U 3n (n = 4, 5) 



ss 



L« = (Ps» - L S1 U 14 - L 62 U 24 - L^U^IAJ, 
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U fiB = C ss -L sl U 15 -L sa U aB -L fi aU 35 -^ 4 U 4S 

Accordingly, the unknowns ao to a 4 can be found from equation (7) below, using the elements !_„,„, U mn and the real 
parts a' §rSl to a' jj>s of the mutual impedance Zq. According to this equation (7), the unknowns ao to a^ can be found 
from Y n obtained from equation (8) given hereinafter and the elements U mn calculated by equation (6). 

a 4 = Y s /U 55 C7) 

where X denotes the sum from n = m+1 to 4, and m = 3, 2, 1,0. Further, U m+1 n denotes the element in the (m+1)th 
row and the n-th column of the upper triangular matrix U t and U m+1 indicates the element in the (m+1)th row and 
the (m+1 )th column of the upper triangular column U. 

Here, Y n is set as shown in equation {6> below. According to equation (8), Y n can be obtained from the elements 
L mn calculated by equation (6) and the values of the real parts a-^ to a'^ determined by the combination of the 
values of i and j. 

m * i|,sm mn n' 

where X denotes the sum from n = 1 torn, and m = 2, 3, 4, 5. Further, L mn represents the element in the m-th row and 
the n-th column of the lower triangular matrix L. 

As is. apparent from the above explanation, the values of the elements L mn and U mn are determined only by the 
values of the five sample frequencies f sl to f^, and are independent of the values of i and j. Accordingly, the values 
of the elements L mn and U mn need only be obtained once at the beginning by solving equation (6), regardless of whether 
the values of the real parts a'ij sl to a- ^ of the mutual impedance Zq change as the values of i and j change. 

The unknowns a 0 to a 4 vary with the values of i and j, and hence with the values of the real parts a'^ sl to a'^ of 
the mutual impedance Zq, but can be obtained by repeating the calculations of equations (7) and (8) according to the 
change of the values of i and j. More specifically, using equation (8), Y-, and Y m are obtained from the elements 
and the values of the real parts a'g s1 to ct' r)AS determined by the combination of the values of i and j, and using equation 
(7), the unknowns ao to a 4 are determined from Y v Y m , and the elements U mn . 

The above explanation also applies for the approximation coefficients b 0 to b 4 for the imaginary parts p'q sl to $'y i s5 
of the approximate equation for the mutual impedance. The same also applies for the mutual admittance and mutual 
reaction. These are apparent from the similarity between the respective approximate equations. 

Figure 4 is a flowchart illustrating the process for calculating the approximation coefficients by the LU decomposition 
method. The series of processing steps S51 to S54 shown in Figure 4 is performed in the approximation coefficient 
calculation processing section 5, and corresponds to step S5 shown in Figure 3 previously explained. The same process 
for the calculation of the approximation coefficients applies, whether the approximate equation is for the mutual im- 
pedance, the mutual admittance, or the mutual reaction. 

As shown in Figure 4, first the elements L,^ and.U mn are obtained from the values of the five sample frequencies 
f s1 to f^ in accordance with equation (6) (S51 ). Next, in accordance with equation (3), specific values of the real parts, 
ot'jj^ to a'q and the imaginary parts, p'q s1 to Pq ^ of the mutual impedance Zq, corresponding to a specific combi- 
nation of the values of i and j, are obtained from the values of the mutual impedance Zq calculated by the sample 
calculation processing section 4 (S52). Then, equations (7) and (8) are calculated (S53), using the values of the ele- 
ments L mn and U mn obtained in S51 and the values of the real parts, a'q s i to cc^sS' and the imaginary parts, to 
P'ijsS' °* the mu tual impedance Zq obtained in S52. In this way, the approximation coefficients ao to a 4 and b 0 to b 4 , 
corresponding to a specific combination of the values of i and j, can be obtained by a relatively simple calculation 
process. The approximation coefficients a 0 to a 4 and b 0 to b 4 thus obtained are for two elements determined by specific 
values of i and j; therefore, the process for obtaining the approximation coefficients a 0 to a 4 and b 0 to b 4 is repeated 
for all combinations of the values of i and j. Then, it is checked whether the approximation coefficients a^ to a 4 and b 0 
to b 4 have been obtained for all combinations of the values of i and j (S54). If the answer is NO, the process is repeated 
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" to obtain the approximation coefficients ao to 84 and b 0 to b 4 for the next combination of the values of . and j. 

In the process of the present invention, therefore, the process step S51 is performed only once, and only the 
processing in steps S52 to S54 is repeated the same number of times as the number of possible combinations of the 
values of i and j. This improves the processing speed. 

Simultaneous Equations Accounting for the Effect of Electromagnetic Field Scattering by Dielectric ■ 

When taking into account the effect of the electromagnetic field scattering caused by the dielectric, in addition to 
the mutual impedance Z° cc between metal segments, the mutual impedances Z° c , d and 7?^ c between metal and 
10 dielectric segments, the mutual impedances Z° d(d and Z^.d between dielectric segments, the mutual admittances 
Y° d d and Yd d d between dielectric segments, the mutual reactions B 0 ^ and B° dc between metal and dielectric seg- 
ments and the mutual reactions B° dd and B* dd between dielectric segments, are obtained. That is, in the actual 
analysis of an electric circuit device, 'all of the mutual impedances, mutual admittances, and mutual reactions are 
obtained Here, the superscript '0" denotes a calculated value in air, the superscript "d" a calculated value in the 
15 dielectric the subscript "C the metal, the subscript "d" the dielectric, the subscript "cC the interrelation from metal to 
metal, the subscript 'd,d' the interrelation from dielectric to dielectric, the subscript "c.d" the interrelation from dielectnc 
to metal and the subscript "d.c" the interrelation from metal to dielectric. 

The'mutual impedance, etc. at sample frequencies are calculated by the sample calculation processing section 4 
using a known calculation process. The mutual impedance, etc. at designated frequencies (e.g., harmonics) within the 
20 specified frequency range are calculated by the approximate calculation processing section 9 in the calculation process 
using the approximate equations of the present invention. Then, the simultaneous equations of the moment method 
are set in the input processing section 1 or the electromagnetic current calculation processing section 10. 

Using the wave source V, present on the printed circuit board, the coefficient l c n denoting the magnitude of the 
electric current flowing in each metal segment, the coefficient l dt0 denoting the magnitude of the equivalent electnc 
25 current flowing on the surface of each dielectric segment, and the coefficient M n denoting the magnitude of the equiv- 
alent magnetic current flowing on the surface of each dielectric segment, the following simultaneous equations of the 
moment method are set, where [] denotes a matrix or vector. 

From the boundary condition that the value of the electric field on the surface of the metal is zero, the following 
simultaneous equations of the moment method are set. 
30 - - 

. [z\ c ][U + [^^^^ 

From the boundary condition that the tangential components of the electric field are equal on both sides of the 
35 boundary surface of the dielectric, the following simultaneous equations of the moment method are set. 

40 From the boundary condition that the tangential components of the magnetic field are equal on both sides of the 

boundary surface of the dielectric, the following simultaneous equations of the moment method are set. 

[B° d . c l[l Cin ! + [B° d , d + B d did ][l d:n ] + [-Y° w - Y d d(d ][M n ] = [0] ^ 

45 

Combining the above, the simultaneous equations of the moment method shown in Figure 5 are set. In setting the 
simultaneous equations of the moment method, if a core realizing a reduction in common mode electric current is 
inserted in the cable, the mutual impedance 7? Ct0 of the cable must be reduced by the impedance ZL of the core. In 
■ that case therefore, the simultaneous equations of the moment method are transformed as shown in Figure 6. 

so Further when there is an earth surface causing reflection, a mirror image of the structure is generated, and elec- 

tromagnetic currents identical in magnitude but opposite in direction to those of the physical structure are assumed for 
the mirror image of the structure. Then, a matrix of the moment method for the physical structure, such as shown in 
Figure 7, is calculated, and also, a matrix of the moment method for the mirror image, such as shown in Figure 8, is 
calculated. Using these two matrix equations, the simultaneous equations of the moment method shown in Figure 9 

55 are set. 

Utilizing the property that the mirror image theoretically has a negative polarity with respect to the physical structure, 
if reflection processing is performed in which an image current is represented by one having a negative polarity with 
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respect to the physical structure current, an increase in the number of unknowns can be avoided, and as a result, the 
calculation time only increases by a factor of about 2. On the other hand, if the image current is treated as an unknown, 
the number of unknowns is doubled, and therefore, the calculation time increases by a factor of about 4. The reflection 
processing, coupled with the approximate calculations, thus achieves a further improvement in the processing time. 

5 The electromagnetic current calculation processing section 1 0 substitutes the mutuaJ impedance, etc. obtained by 

approximate calculations into the above-given simultaneous equations of the moment method, and by solving the 
simultaneous equations, calculates the electric current flowing in the metal on the printed circuit board, the equivalent 
electric current and equivalent magnetic current flowing on the surface of the dielectric on the printed circuit board, the 
electric current (including the common mode electric current) flowing to the cable, wire, and leads, and the common 

io mode electric current flowing to the metal housing and metal structure. 

Accuracy of Approximate Calculations 

To verify the validity of the present invention, the mutual impedance values obtained by the approximate calcula- 
15 tions of the present invention were compared with those obtained without using approximate calculations, and the 
errors were examined. 

Figure 10 shows sample frequencies and also a frequency for which the impedance was calculated by approxi- 
mation. In Figure 10, values of wavenumber k corresponding to the respective frequencies are also shown. 

Five sample frequencies were taken, so that n = 5. Therefore, the approximate equation, Zjj = exp(-jkr)[ao + a-,f 2 
20 + agf 4 + a 3 f 6 +.a 4 f B + j(b(/f + b.,f + b 2 f 3 + b 3 f 5 + b 4 f 7 )], was used for the calculation of the mutual impedance 2^. The 
five sample frequencies, f s1 , f^, f^, f^, and f^, were chosen to be 30 MHz (3e + 007 or 3 X 107), 250 MHz, 500 MHz, 
750 MHz, and 1 GHz. This setting is for the frequency range of 30 MHz to 1 GHz. The mutual impedances at the 
five sample frequencies are obtained by the sample calculation processing section 4 using a known calculation process 
used in the moment method. 

25 The frequency f for which the mutual impedance 2^ was to be obtained by approximate calculation was set at 1 00 

MHz (le-+ 008 or 1 x 10 8 ). In this comparison test, the monopole structure shown in Figure 11 was assumed. The 
comparison test assumed that the monopoles (1)~and (2) shown at the extreme left in the figure were fixed. On the 
other hand, monopoles (3) and (4) were chosen to be one of #1 to #6 in the figure, the monopole-to-monopole distance 
r thus being variable. In Figure 11 , it is shown as if the monopoles were placed at equally spaced intervals; actually, 

30 the distance from the monopoles (1) and (2} to the dipole #1 is 0.0001 m (meter), the distance to the dipole #2 is 0.001 
m, the distance to the dipole #3 is 0.01 m, the distance to the dipole #4 is 0.1 m, the distance to the dipole #5 is 1 m, 
and the distance to the dipole #6 is 10 m. Wire diameter, i.e., monopole diameter, was 2um.The length of each monopole 
was 0.03 m. 

Figure 1 2 shows calculated values of the mutual impedance 2^ when the dipole-to-dipole distance r = 0 m. Since 
35 the distance r = 0 m, Figure 1 2 actually shows the values obtained by calculating the self-impedance of the monopoles 
(1 ) and (2) using the approximate equation. In this patent specification, the concept of mutual impedance includes self- 
impedance. 

In the approximation coefficient section of Figure 12, c 0 , c lP c 2 , c 3 , and c 4 respectively indicate the calculated 
values of the approximation coefficients aQ, a 1( a 2 , a 3 , and a 4 for the real part of the above approximate equation. 
40 Likewise, m 0l m 1 , m 2 , m 3 , and m 4 respectively indicate the calculated values of the approximation coefficients b 0 . b-, , 
b 2 , b 3 , and b 4 for the imaginary part of the same approximate equation. In Figure 12, Cq, a,, c 2 , c 3 , and c 4 are different, 
and therefore to be distinguished, from the approximation coefficients used in the earlier given approximate equation 
for the mutual reaction. 

The approximation coefficients c 0 to c 4 and m 0 to m 4 are obtained in the approximation coefficient calculation 
4S processing section 5 by using the five sampling frequencies, etc. 

Figure 1 2 also shows values each obtained by multiplying each approximation coefficient by a power of wavenum- 
ber k. For example, c,*k A 2 indicates the value obtained by multiplying the approximation coefficient c, by the second 
power of wavenumber k. The value of wavenumber k was obtained by substituting the frequency f = 100 MHz to be 
analyzed into the equation k = 2kVc. Multiplying an approximation coefficient by a power of wavenumber k has the 
. 50 same meaning as multiplying the approximation coefficient by a power of the frequency f. The same also applies in 
Figures 13 to 18. 

The exponent part section shows the value of erfa in the approximate equation. In the exponent part section, the 
value of the real part and the value of the imaginary part are shown in parentheses, the former before the comma and 
the latter after the comma. The approximate value section shows the value of the mutual impedance calculated by the 
55 approximate equation of the present invention. In the approximate value section, the value of the real part and the 
value of the imaginary part are shown in parentheses, the former before the comma and the latter after the comma. 
The above values are obtained in the approximate calculation processing section 9. The same also applies in Figures 
13 to 18. 
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The moment method section shows the value of the mutual impedance calculated without using the approximate 
equation In the moment method section, the value of the real part and the value of the imaginary part are shown in 
parentheses, the former before the comma and the latter after the comma. These values were obtained using the prior 
art electromagnetic field intensity calculation apparatus. The same also applies in Figures 1 3 to 1 8. 

The error section shows the error between the mutual impedance value obtained without using the approximate 
equation and the mutual impedance value calculated by the approximate equation of the present invention. In the error 
section, the error of the real part and the error of the imaginary part are shown in parentheses/the former before the 
comma and the latter after the comma. . 

As can be seen from the comparison of the values shown in the moment method section and the approximate 
value section, the values are very nearly equal between the two sections. As shown in the error section, when the error 
is calculated to five decimal places of significant figures, the error becomes 0%. This shows an extremely good ap- 
proximation. . « ««- j 
- Figure 1 3 shows the values when the dipole-to-dipole distance r = 0.0001 m, Figure 1 4 when r = 0.001 m, and 
Fiqure 1 5 when r = 0 01 m As shown, in these cases also, the error, when calculated to five decimal places of significant 
figures, is 0%. It is therefore safe to say that in the case of a frequency of 100 MHz, there are no errors when the 
distance is 0.01 m or less. , , r 1D 
Figure 16 shows the values when the dipole-to-dipole distance r = 0.1 m, Figure 17 when r = 1 m, and Figure 18 
when r = 1 0 m In the cases shown in Figures 1 6 to 1 B, errors occur as indicated. In the case of the frequency 1 00 
MHz, the error peaks at r = ai m. As can be seen from the figures, the error does not increase wrth increasing dipole - 
distance r, but converges. This also shows that an extremely good approximation is achieved. 

Proof that Mutual Impedance, etc. can be Approximate d bv Polynomials in Frequency 

We will now prove that the mutual impedance can be approximated using the approximate equation. For an ex- 
planation about the approximate equation tor the mutual impedance, we consider the monopoles shown in Figure 1 9. 
In Figure 19 thick lines show the monopoles and dotted lines show the shapes of expansion functions. Monopoles 
(1 ) (2) and monopoles (3), (4) are respectively located in two parallel planes separated by a distance h, and the angle 
between them is denoted by <ft. More generally, the tilt angle between monopole (1) and monopole (3) is denoted by 
A, and the tilt angle between monopole (1 ) and monopole (4) by <t> 2 .- 

The mutual impedance Z between the two monopoles whose expansion functions are J, and J 2 is expressed by 
the following general equation, where © is the angular frequency, r is the distance between the monopoles, Pl = -1/jco 
X 3 /at and p 2 = -1 /jo x aJa^t- 1 n the equation below, th e integration is performed with respect to s so that the equation 
can be applied as a general equation not only when the monopole shape is linear {a wire) but when it is a surface (a 
surface patch). 

• z = 3« W^r Ji j2COS * f~ + /4^r PiPz ~ — ]ds 

The expansion functions J, and J 2 ' differ according to the type of the moment method, but the following proof is 
applicable regardless of the type of the moment method. More specifically, the expansion, functions may be those for 
a sinusoidal current (in the case of piecewise sinusoidal moment method), a triangular current, or a pulse function 
current distribution Therefore, the mutual impedance calculation process by the approximate equation of the present 
invention can be applied regardless of the type of the moment method. The following proof is also applicable whether 
the monopole shape is of a surface (surface, patch) or linear (wire). Accordingly, the mutual impedance calculation 
process by the approximate equation of the present invention can be applied regardless of the monopole shape. 

. We will now prove that the mutual impedance can be approximated by a polynomial with terms each involving 
raising the frequency f to an n-th power (n = -1, 0, 1, 2, ...). In the following proof, however, a more general form, 
wavelength k, may be used as necessary. Here, k = <h/c, where c is the velocity of light. 

The following prool specifically concerns the piecewise sinusoidal moment method, but rt will be appreciated that 
the same proof is also applicable for other modes of the moment method that use a triangular current or a pulse function 
current distribution as the expansion function. 

In the piecewise sinusoidal moment method, the expansion functions for the monopoles (1 ) to (4) of Figure 1 9 can 

be expressed as : 

Electric Monopole (1) J, = sink(z-z 0 )/sinkd-, 
Electric Monopole (2) J n = sink(z 2 -z)/sinkd 2 
Electric Monopole (3) J 2 = sink(t-t 0 )/sinkd 3 
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Electric Monopole (4) J2 = sink(t 2 -t)/sinkd 4 

where d-, is the length of the monopole (1 ), is the length of the monopole (2), d 3 is the length of the monopole (3), 
and d 4 is the length of the monopole (4). 
s Using these expansion functions, the mutual impedance Z 13 between the monopoles (1>and (3>and the mutual 

impedance Z 14 , between the monopoles (1) and (4) are obtained first. The mutual impedances Z, 3 and Z 14 are given by 

*Z 13 = -2 . • ^ J J[sink(z - z 0 )sink(t - tojcos*! 

t 0 z 0 e" jkr 
- coskjz - z 0 )cosk(t - t 0 )] — dzdt 

t 2 z i 

15 2 14 = ■ U o / J[sink(z - z 0 )sink(-t + t 2 )cos<J> 2 

14 4nsinkd!Sinkd 4 L 1 

+ cosk(z - 2 0 )cosk(-t + t 2 ) ] — dzdt 

20 The coefficients for the mutual impedances Z 13 and Z 14 in the above equations are expressed as 

\m = k =\ a — 

47isinkd 1 sinkd 3 J 4it sinkd 1 sinkd 3 1 A 1 

25 

lm = jCy k =\ a — 

47isinkd 1 sinkd 4 J 4n sinkd 1 sinkd 4 A 2 

30 where a - cu/4rc 

Further, the monopole-to-monopole distance r can be approximated, using r 0 = (Z-, 2 + t, 2 - 2z 1 t 1 cos4> + h 2 ) 1 ' 2 , by 

,2 ,2 2..1/2 .2 A ,1/2 .„ . . 2.1/2 
r = [r 0 + ( r - r o)l =( r o +A > =r 0 ( 1+A/r o) 



35 



^r 0 (1 + A/2r 0 2 -A 2 /£r 0 4 + .„) 
= r 0+ A/2r 0 -A 2 /8r 0 3 + ... = r 0 + d 



40 Hence the mutual impedances Z 13 and Z 14 are expressed as 



55 



k tlZl 
Z u = jet — — exp(-jkr 0 )| /[sink(z - z 0 )sink(t - tojeos*! 



4S , t 0 z 0 



e -J" 



- cosk(z - z 0 )cosk(t - t 0 )] — dzdt 



so v t 2 z i 

= D a -vt — ex P("j kr o)/ /[sink(z - z 0 )sink(-t + t 2 )cos(t> 2 

tlZ ° e -J" 
+ cosk(z - z 0 )cosk(-t + t 2 ) ] dzdt 



where = sinkd^inkd^j 

Then, setting z - z Q = u, t - ^ = v,. and w = -t + t 2 . the equations are put in simplified form. Further, assuming that 
the monopole length is short, the following polynomial approximation is performed. 



16 



EP 0 778 533 A2 

sinku = ku - (ku) 3 /6, cosku = 1 - (ku) 2 /2 

Calculation example: 

u = 0.1A - sin(0.2n) = 0.587785 

^ ku - 0.6283 (+6.89%) 

^. ku - (ku) 3 /6 = 0.586958 (-0.14%) 

cos(0.2n) = 0.809016 

= 1 - (ku) 2 /2 = 0.802608 (-0,79%) 

u= 0.12X - sin(0.24rt) =0.684547 

,ku = 0.75398 (+10.14%) . 

' \ ku - (ku)V6 = 0.682542 (-0.29%) 
cos(0.24u) = 0.728969 
^ 1 -(ku) Z /2 = 0.715755 (-1.81%) 

where X is the wavelength, or 



From the above, the real part of the mutual impedance Z 13 is expressed as lollows. In the equation below, exp( 
jkr 0 ) is omitted. 
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k ttZi 

a — J J[s inku sinkv cos^t 

Ai toZo 

sinkd 

-cosku coskv} dzdt 

r 

k hz, <ku> 3 (kv) 3 

a — S f CCku 3Ckv Dcostf , 

A. t,2o 6 6 

. (ku) 1 (kv)\ (kd>! 1 

-CI DO 3 3 Ckd ] -d2dt 

2 2 • 6 r 



k t.z, uv 3 + u 3 v u 3 v 

a — J J" (Cuvk^ k< + — k s :co S9 5 

A, toZo • 6 ' 36 



uHv 1 uV ' (kd) 3 ' 1 : 

_q. k *x 4<<] ) Ckd- — — D — dzdl 

2 4 -6 r 



1 l,2i • uv 3 + u 3 v uvd 3 

a - J / CCuvdkM d+ — )k' 

Ai Uzo • 6-- 6 



uV (uvH u 3 v>d 3 . u 3 v 3 d 3 
+ ( d+ )k 8 k ,0 ]cos 4>\ 

- 36 • 36 216 . 



uHv 1 d' uV "* v * ,, t 

-dk'+( d+ — )k<-( d+ d 3 )k 

2 6 4 12 



uV 1 • 

+ dV ] ) d2dt 

24 r 
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Replacing the coefficients for the terms with k* etc. by P, to P 8 to simplify the equation in form, the above equation 
can be written as 

= [oJA,][(P, k 4 - P 2 k 6 + P 3 k 8 - P 4 k 10 )cos<l> 1 - P 5 k 2 + P 6 k 4 - P 7 k 6 + P 8 k 8 ] 

Further since ^ is a sin function, this can be approximated as A, = sinkd n sinkd 3 = d^k*. Substituting this ap- 
proximate value of A, and inserting the previously omitted exp(-jkr 0 ) back into the equation, the above equation can 
be written as 

' R, = [aexp(-jkr 0 y{d l d 3 k 2 )][(P 1 k 4 - P 2 k 6 + P 3 k 6 P 4 k 1 Vs<h - P 5 k 2 + P 6 k 4 - P 7 k 6 + P 8 k 8 ] 
. = [aexpHkroVCd.dgHKP^ 2 - P 2 k 4 + P 3 k 6 - P 4 k 8 )cos<t> 1 - P s + P 6 k 2 - P 7 k 4 + P 8 k 6 ] 

The real part R n of the mutual impedance Z 13 has thus been obtained. 

Next, the imaginary part I, of the mutual impedance Z 13 is expressed as follows. In the equation below, exp(-]kr 0 ) 

is omitted. 

* k tizt 

[ ,= a — S S [sinku sinkw cos <p i 

Ai toZo 

coskd 

4cosku coskw] dzdt 



30* 



- a 



35 At toZo 



k t.zi (ku) 3 Ckw) 3 

J* J.CCku ][kw ]cos^i 



40 



45 



.50 



55 
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(ku) J Oro) 1 (kd) J 1 ' 

[1 3D ] 3 [1 ] — didt 

2 2~ 2 r 



k tiZl uw 3 + u 3 w u 3 w 3 

a — S S C Cuwk 7 kH — k*3cos <6 . 

Ai to2o 6 35 



u'+w 3 u'w 1 (kd> 2 I- 

_ 1+ k * k 4 ] 3U 33— dzdt 

2 4 2 r • 



c t,z, uw 3 i Vw uwd 1 
— / / [ Cuwk 3 -( ■ + — >k s 

Al taZo 6 2 

uV uff'f u 3 w u 3 w 3 

+ ( + d 3 )k 7 d 3 k'3cos0 i 

36 12 72 . 



uHv? 3 d J • u 2 w 3 uHw 3 

-kf( + )k 3 -( + d 3 )k s 

2 2 4 4 



u V 1 

+ — dV 3 dzdt 

8 r 

Replacing the coefficients for the terms with k 3 , etc. by Q, to Q a to simplify the equation in form, the above equation 
can be written as 

' It = [a/A, ][(Q 1 k 3 - Q 2 k 5 + Q 3 k 7 - Q^Jcosc^ - Q s k + Q 6 k 3 - Q 7 k 5 + Q 8 k 7 ] 

Further, substituting the approximate value of A, = sinkc^sinkda = d n d 3 k 2 and inserting the previously omitted exp 
(-jkr 0 ) back into the equation, the above equation can be written as 

I, = [aexp(-ikr 0 )/(d 1 d 3 k 2 )][(Q 1 k 3 - Q 2 k 5 + Q 3 k 7 - 
Q 4 k 9 )cos* 1 - Q s k + Q 6 k 3 - Q 7 k 5 + Q 8 k 7 ] = ^exp^r^d^)]^ k - Q 2 k 3 + Q 3 k 5 - Q 4 k 7 )cos(> 1 - Q 5 /k 

+ Q 6 k-Q 7 k 3 +Q a k 5 ] 
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The imaginary part I, ol the mutual impedance Z, 3 has thus been obtained. 
Hence, the mutual impedance Z-, 3 is expressed as 

Z 13 = R 1 -jl, =[ae X p(-ikr 0 y(d 1 d3)][(P 1 k 2 -P 2 k 4 + P 3 k 6 - P 4 k 8 )cos*, - P s + P 6 k 2 - P 7 k 4 + P s k 6 ] + 
KaexpHkr^td^I^k - Q 2 k 3 + Q 3 k 5 - a/jcos^ - Q 5 /k - Q 6 k - Q 7 k 3 + Q 8 K S ] 

The mutual impedance Z, 4 . can also be obtained in like manner. The equation expressing the mutual impedance 
Zl4 is lundamentally the same as that for the mutual impedance Z n3 , except that the terms not multiphed by cos* are 
reversed in sign. That is. the mutual impedance Z 14 , is expressed as shown below. Here, R and S are used as the 

coefficients corresponding to P- and Q. 

Z 14 = [aexpHkr^d,)]^ k 2 - R 2 k 4 + R 3 k S - R 4 k 8 )cos* 2 + * 5 - + *7»< 4 ' ^ + 
HaexpHkroW^dJKS, k - S 2 k 3 + S 3 k 5 - S 4 k 7 )cos$ 2 + S 5 /k - S,k + S 7 k 3 -S 0 k 5 ] 

Hence, the mutual impedances (Z 13 +.Z 14 ) are generally expressed as 

Z 13 + Z 14 = exp(-jkr 0 )[ (C 0 + C, k 2 + C 2 k 4 + C 3 k 6 + C 4 k 8 ..... j^k 1 + C 6 k + C ? k 3 + C s k 5 + C 9 k 7 + ..)] 

The mutual impedances (Z 23 + Z, 4 ) also can be obtained in like manner. Z 23 is the mutual impedance between 
the monopoles (2) and (3) and Z 24 the mutual impedance between the monopoles (2) and (4). 

According*! the mutual impedances (Z 13 + Z 14 + Z^ + Z 24 ) can be expressed by a po.ynom.al ,n wavenumber k 
(i e " frequency f),- similarly to the above equation. . ftrnnttan „, 

' For the mutual admittance, it is possible to prove that the approximation can be done by a polynom.a! of frequency, 
in a manner similarto the mutual impedance, and the approximate equation is identical tothatforthe mutual impedance, 
as previously described. Accordingly, the proof is not given here. ranamta „ a tinn 

It will be proved that the mutual reaction can be approximated using the approximate equat.on. For an explanation 
of the approximate equation for the mutual reaction, the monopoles shown in Figure 1 9, where, of the monopoles (1 ) 
to (4), the monopoles (3) and (4) are regarded as magnetic current monopoles, will be considered. 

First, a magnetic field produced by an electric current source is considered. . 
" The expansion functions for the monopoles (1) to (4) can be expressed as 

Electric current monopole (1 ) = sinkfz-ZoVsinkd, 
Electric current monopole (2) J 2 = sink(z 2 -z)/sinkd 2 
Magnetic current monopole (3) = sink(t-to)/sinkd 3 
Magnetic current monopole (4) M 4 = sink(t 2 -t)/sinkd 4 

Further, unit vectors can be expressed as shown below. 
Unit vector: 

Monopole (1)(s X) Sy,s 2 )= (0,0,1) . 
Monopole (3) (V t y , t z ) = (sintt),, 0, cos^) 
Monopole (4)- (sin$ 2 , 0, cos4> 2 ) 
p direction (v x , v y , v 2 ) = (x/p, h/p, 0) 
<t> direction (w X! w y , w 2 ) = (-h/p, x/p t 0) 

When the electric current source exists only on the z-coordinate, the electromagnetic field becomes circularly sym 
metrical (independent of the ^-coordinate), and only H0 exists, which can be expressed as 



u, [i dZ dp 



21 



10 



1S 



20 



25 



30 



Ap = 0, hence 
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u. dp 



Next, we consider the reaction. 

In Figure 1 9, since the tangential components of the magnetic field on the monopoles (3) and (4) are given by Ht 
= -h/p x H<{>sin<> t the reaction B can be expressed as - 

B = -h J M ^-sintt) hdt 
t P 

= J}_J M _J_ J>± sincb dt 
\x J t p dp 

Here, the following conditions are substituted into the above equation. 



-ikr 

z 4n J r 



- 2 .2.1/2 
p = {X + h ) 



2 2 1/2 2 2 2 1/2 

r = {p + {z - z') } = (z + 1 - 2ztcos § + h ) z': z coordinate of monopble (3), (4) 



dr p 
dp ~ r 



35 dA z _ u. f T c e 



a P 4— .J J -a7 -T- dz 



= 7 J <-J* -V 31tr - -e T ^- ,k *)d» 



z r 



Hence, the reaction B can be expressed as 



B = / J J M — -f— -^-sin(}> dz dt 

4n J o od . r 



45 D 4rc J J p c>p , 



4 " t z 

50 Further, since r = r 0 + d, the reaction B can be expressed as 



B * hsin» o - Jfc ro f s M j _ik L-Je' 3 " dz dt 

4n t z r r 

55 

Accordingly, the reaction B 13 between the monopoles (1 ) and (3) and the reaction B 14 between the monopoles (1 ) and 
(4) can be expressed as 
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B, 3 - 



hsin< * > i e" ikT ° f f sink(z - z 0 ) 



13 4nsinkd! sinkd 3 *. z 

sink(t - t 0 ) I- -^r p-l e " 3 " dz dt 

B = hsinfrz e - JMO / f sink(z - z 0 ) 

Bl * 4rcsinkd, Sinkd t t z > 

sink(t 2 .- t) [- - ^-]e" J " dz dt 

Then, setting z - z$= u, t - 1 0 = v, and w = -t + 1 2 , the equations are put in simplified form. Further, assuming that the 
monopole length is short, the following polynomial approximation is performed. That is, sinku*= ku - {ku) 3 /6 and cosku 
= 1 ■- (ku) 2 /2 

From the above, the second term on the right-hand side of the mutual reaction B, 3 equation can be written as 

1 ' „ 
/ /sinku sinkv [ - — — ]e dz dt 

t z r 3 

(ku) 3 ■ (kv) 3 

=-J J(ku ][kv-- J 

t z 6 6 

1 • ■ 

[coskd- j sinkd] dz dt 

r 

(ku) 3 - r (kv) 3 (kd) 2 ; 

=-//[ [ku- — ][kv -— — 3 

t z 6 . 6-2 

(ku) 3 (kv) 3 ' 
- j . Cku -.- — ][kv - ] 
6 6 

(kd) 3 ' 1 
[kd - — ] ]— dzdt 
6 r 3 
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u 3 y+uv 3 + 3uvd 2 
■f Si uvk 2 k< ■ 

t 2 6 

u 3 vWvdH3uv 3 d 2 u 3 v 3 d 2 

+ : : k 6 -- k 8 

36 72 



u 3 vd+uv 3 d+uvd 3 
. -j [uvdk 3 k s 



uVd+u 3 vd 3 +uv 3 d 3 uVd 3 1 . 

+ = k 7 k 9 ) ] — 62. dt 

36 216 r 3 

Replacing the coefficients for the terms with k 2 , etc. by P, to P 8 to simplify the equation in form, the above equation 
can be written as 

P A 2 + P2 k4 + p 3 k6 + p 4 kB + i( p 5 k3 + p & k5 + p 7 k? + Pek 9 ) Likewise, the first term on the right-hand side of the 
reaction &,3 equation can also be written as 

jk 

J Jsinku sinkv ( ]e dz dt 

t 2 r 2 

u 3 v+uv 3 + 3uvd 2 

=-//'[ uvk 2 : k 4 

t 2 6 • 



uV*3u 3 vd 2 +3uv 3 d 2 s u 3 v 3 d 2 ^ 
36 - 72. 
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u 3 vdTUV 3 d+uvd 3 
-j (uvdk 3 - : k 5 



u 3 v 3 d+u 3 vd 3 +uv 3 d 3 , u 3 v 3 d 3 qi ^ jk 

+ k 7 k 9 ] ] — d2 dt. 

36 216 r 2 



u 3 vdTuv 3 d+uvd 3 

=- 1 S ( uvdk 4 k s 

t z - 6 



• uVd+u 3 vd 3 -hjv 3 d 3 u 3 v 3 d 3 

+ : k a k 10 

36 216 



u 3 v+uv 3 + 3uv.d 2 
+j[uvk 3 - — : k s 



uV+3u 3 vdH3uv 3 d 2 uVd 2 1 

+ — k 7 k 9 ] ]— - dz dt . 

36 72. r 2 

Replacing the coefficients for the terms with k 2 , etc. by'Q, to Q 8 to simplify the equation in form, the above equation 
can be written as ' . - 

Qy + Q 2 k 6 +Q 3 k 8 + Q 4 k 10 + J(Q 5 k 3 +Q 6 k 5 + Q 7 k 7 + Q 8 k 9 ) 

Accordingly, the mutual reaction B 13 is expressed as 

B 13 = (h t sin<t>/4iisinkd 1 sinkd 3 )e*' k, ° X [P^ 2 + P 2 k 4 + P 3 k 6 
+ P 4 k 8 + j(P 5 k 3 + P 6 k 5 + P 7 k 7 + P 8 k 9 ) + Q n k 4 + Q 2 k 6 + Q 3 k 8 + Q 4 k 10 + j(Q 5 k 3 + Q 6 k 5 + Q 7 k 7 + Q 8 k 9 )] 

Substituting the approximate value, sinkd 1 sinkd 3 H d 1 d 3 k 2 J into the above equation, the mutual reaction B 13 can be 
expressed as 
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B 13 = (hsin<t> 1 /4nd 1 d 3 )e~ pa0 x [P 1 + (P 2 + Q^k 2 + (P 3 + 
Q 2 )k 4 + (P 4 + Q 3 )k 6 + Q 4 k 8 + j[(P 5 + Q 5 )k + (P 6 + Q 6 )k 3 + (P 7 +' Q 7 )k 5 + (P 8 + Q 8 )k 7 ]] 
= e" ito0 [R 1 + Rjk 2 + R 3 k 4 + R 4 k 6 + R 5 k 8 + j(R 6 k + R 7 k 3 + R 8 k 5 + \%k 7 )] 

The mutual reaction B 14 , can be obtained in (ike manner. That is, the mutual reaction B 14 , is expressed as 

B 14 = e" ila0 [S 1 + S 2 k 2 + S 3 k 4 + S 4 k $ + S s k 8 + ](S 6 k + Syk 3 + S 8 k 5 + S g k 7 )] 
Accordingly, the mutual reactions (B^ + B 14 ) can be written as 

B 13 + B 14 = e'^fC, + C 2 k 2 + C 3 k 4 + C 4 k S + C 5 k 8 + j(C 6 k + C 7 k 3 + C B k s + C g k 7 }] 

The mutual reactions (Bgg + B 24 ) also can be obtained in like manner. B^ is the mutual reaction between the 
monopoles (2) and (3) and Bg 4 the mutual reaction between the monopoles (2) and (4). 

Accordingly, the mutual reactions (B^ + B 14 + B^ + B24) can be expressed by a polynomial in waven umber k (i. 
e., frequency f), similarly to the above equation. 

Extension ot the Approximation Coefficient 

Calculation Process 

As previously described, the process in the approximation coefficient calculation section 2 has been such that the 
mutual impedance, etc. are calculated using known methods in the sample calculation processing section 4 without 
using approximate equations and, using the calculated results, simultaneous equations are solved in the approximation 
coefficient calculation processing section 5 to obtain the approximation coefficients. However, the process for obtaining 
the approximation coefficients in the approximation coefficient calculation section 2 is not limited to the illustrated 
process. 

For example, the approximation coefficients may be computed by direct calculations. As indicated in the proof 
given for the deriving of the approximate equation for the mutual impedance, the approximation coefficients are a 
function of distance. Briefly described, the coefficients C 0 , etc. at the stage immediately before deriving the approxi- 
mation coefficients are expressed by the coefficients, P v etc., Q.,, etc., and S v etc. at the immediately preceding stage. 
These coefficients P 1t etc. are in turn expressed by u, v, w, d, etc. Since the values of u, etc. define the monopole-to- 
monopole distance r, these values represent the distance; since the monopole-to-monopole distance r can be defined 
as described above, its value can also be obtained. This shows that the approximation coefficients can be obtained 
by direct calculations. 

This also holds true for the mutual admittance and mutual reaction. Accordingly, the approximation coefficient 
calculation section 2 may be configured so that the plurality of approximation coefficients used in the approximate 
equations expressing the mutual impedance, mutual admittance, and mutual reaction can be computed by direct cal- 
culations. 



Claims 

1. A method of calculating an electromagnetic field intensity, comprising the steps of: 

a) calculating approximation coefficients for each of approximate equations that respectively yield approximate 
values of electromagnetic characteristic values, between a plurality of elements constituting an electric circuit 
device and of the elements themselves, at an arbitrary frequency; and 

b) calculating the approximate values of the electromagnetic characteristic values between the elements and 
of the elements themselves at a designated frequency in accordance with the approximate equations having 
the approximation coefficients calculated in step a), thereby making it possible to evaluate the electromagnetic 
field intensity for the designated frequency. 
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2. A method according to claim 1 , further comprising the steps of: 

c) calculating an electric current flowing on each element by solving simultaneous equations obtained by sub- • 
stituting the approximate values of the electromagnetic characteristic values calculated in step b) into equations 
respectively describing boundary conditions for the elements; and 

d) calculating the electromagnetic field intensity for the designated frequency from the electnc currents re- 
spectively flowing on the elements calculated in step c). 

3. A method according to claim 1 , wherein the step a) includes the substeps of: 

i) calculating the electromagnetic characteristic values between the elements and of the elements themselves 
from geometric data of the elements at each of a plurality of sample frequencies; and 

ii) calculating the approximation coefficients by solving simultaneous equations obtained by substituting the 
Electromagnetic characteristic values calculated at the plurality of sample frequencies into each of the approx- 
imate equations. 

4. A method according to claim 3. wherein the substep a) ii) includes the substeps of: 

calculating matrix elements in a lower triangular matrix and an upper triangular matrix from the sample fre- 
' quencies said lower triangular matrix and upper triangular matrix being derived from a matrix containing pow-. 
ers of the sample frequencies as matrix elements by being decomposed as a product of a lower tr.angutar 
matrix, with diagonal elements of 1, and an upper triangular matrix; and 

calculating the approximation coefficient from the electromagnetic character.st.c values and the matrix ele- 
ments in the lower triangular matrix and the upper triangular matrix. 

5 A method according to claim 1 . wherein the electromagnetic characteristic values include mutual impedances 
' " between the elements constituting the electric circuit device and of the elements themselves. 

6. A method according to claim 5, wherein 

the elements constituting the electric circuit device include conductor elements and dielectric elements, and 
the electromagnetic characteristic values further include mutual admittances between the dielectric elements 
and of the dielectric elements themselves, mutual reactions between the conductor and dielectric elements, 
and mutual reactions between the dielectric elements and of the dielectric elements themselves. 

7. a method according to claim 6, further comprising the steps of: 

c) calculating an electric current, an equivalent electric current, and an equivalent magnetic current flowing in 
each element, by solving simultaneous equations obtained by substituting the approximate values of the elec- 
tromagnetic characteristic values calculated in step b) into equations respectively describing boundary con- 
ditions for the elements; and 

d) calculating the electromagnetic field intensity for the designated frequency from the electnc currents, the 
equivalent electric currents, and the equivalent magnetic currents respectively flowing on the elements, cal- 
culated in step c). 

B A method according to claim 5, wherein the approximate equations for the mutual impedances Zj, are 2 r , = exp(- 
Ma, + + a 2 f< I a 3 fe + a 4 f* + ... + KM + b,f + b 2 f3 + b 3 f* + b 4 f + ...)]. where is the frequency, k is wavenumber 
where k = 2irf/e and c is the velocity of light, r is an element-to-element distance, ] is an imaginary unit, and a* a, 
and b 0 , b v ... are the approximation coefficients. 

9 A method according to claim 6, wherein the approximate equations for the mutual admittances Yg are Y„ = exp(- 

where k = 2irf/c and c is the velocity of light, r is an element-to-element distance, j is an imaginary unit, and a* a, 
and b 0 , b v ... are the approximation coefficients. 

10 A method according to claim 6, wherein the approximate equations for the mutual reactions Bq are = exp(-jkr) 
[cn + cJ* + c 2 f* + 4* + c 4 fB + ... + i(d 0 f + d,f3 + d 2 f* + d 3 f? + d 4 f9 + ...)]. where t is the frequency, k « wavenumber 
where k = 2nf/c and c is the velocity of light, r is an element-to-element distance, j is an imaginary unit, and C* 
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c 1t ... and d 0 , d v ... are the approximation coefficients. 



11. A method according to claim 3, wherein the number of sample frequencies is equal to the number of terms in a 
real part and in an imaginary part of each of the approximate equations. 

5 

12. A method according to claim 11 , further comprising the step of determining values of the plurality of sample fre- 
quencies based on a specified frequency range and on the number of sample frequencies. 

13. An apparatus for calculating an electromagnetic field intensity, comprising: 

10 

means for calculating approximation coefficients for each of approximate equations that respectively yield 
approximate values of electromagnetic characteristic values, between a plurality of elements constituting an 
electric circuit device and of the elements themselves, at an arbitrary frequency; and 

means for calculating the approximate values of the electromagnetic characteristic values between the ele- 
js ments and of the elements themselves at a designated frequency in accordance with the approximate equa- 

tions having the approximation coefficients calculated by the approximation coefficient calculating means, 
thereby making it possible to evaluate the electromagnetic field intensity for the designated frequency. 

14. An apparatus according to claim 1 3, further comprising: 

20 

means for calculating an electric current flowing on each element by solving simultaneous equations obtained 
by substituting the approximate values of the electromagnetic characteristic values calculated by the approx- 
imate value calculating means into equations respectively describing boundary conditions for the elements; 
and 

25 means for calculating the electromagnetic field intensity for the designated frequency from the electric currents 

respectively flowing on the elements, calculated by the electric current calculating means. 



15. An apparatus according to claim 13, wherein the approximation coefficient calculating means includes: 

30 means for calculating the electromagnetic characteristic values between the elements and of the elements 

themselves from geometric data of the elements at each of a plurality of sample frequencies; and 
means for calculating the approximation coefficients by solving simultaneous equations obtained by substi- 
tuting the electromagnetic characteristic values calculated at the plurality of sample frequencies into each of 
the approximate equations. 

35 

16. An apparatus according to claim 15, wherein the approximation coefficient calculating means includes: 

means for calculating matrix elements in a lower triangular matrix and an upper triangular matrix from the 
sample frequencies, said lower triangular matrix and upper triangular matrix being derived from a matrix con- 
40 taining powers of the sample frequencies as matrix elements by being decomposed as a product of a lower 

triangular matrix, with diagonal elements of 1, and an upper triangular matrix; and 

means for calculating the approximation coefficient from the electromagnetic characteristic values and the 
matrix elements in the lower triangular matrix and the upper triangular matrix. 



<s 17. An apparatus according to claim 13, wherein the electromagnetic characteristic values include mutual impedances 
between the elements constituting the electric circuit device and of the elements themselves. 

18. An apparatus according to claim 17, wherein 

so the elements constituting the electric circuit device include conductor elements and dielectric elements, and 

the electromagnetic characteristic values further include mutual admittances between the dielectric elements 
and of the dielectric elements themselves, mutual reactions between the conductor and dielectric elements, 
and mutual reactions between the dielectric elements and of the dielectric elements themselves. 

55 19. An apparatus according to claim 18, further comprising: 

means for calculating an electric current, an equivalent electric current, and an equivalent magnetic current 
flowing on each element, by solving simultaneous equations obtained by substituting the approximate values 



28 



EP 0 778 533 A2 

of the electromagnetic characteristic values calculated by the approximate value calculating means into equa- 
. ■ tions respectively describing boundary conditions tor the elements; and 

means for calculating the electromagnetic field intensity for the designated frequency from the electric currents, 
the equivalent electric currents, and the equivalent magnetic currents respectively flowing on the elements 
thus calculated. 

20 An apparatus according to claim 17, wherein the approximate equations for the mutual impedances 2^ are = 
exp(-jkr)[ao + a^ + a 2 f^ -+ a 3 f e + a 4 f* + ... + j(b</f + b n f + b 2 P +■ b 3 f* + b 4 f7 + ...)], where f is the frequency, k is 
wavenumber where k = 2tcf/c and c is the velocity of light, r is an element-to-element distance, j is an imaginary 
unit, and ao, a 1t ... and b 0 , b v ... are the approximation coefficients. 

21 An apparatus according to claim 18. wherein the approximate equations for the mutual admittances Y Vi are Y fj = 
exp(-jkr)[ao + aj* + a 2 H + a 3 f 6 + a 4 f 8 + ... + j(b</r + b,f + b 2 f3 + b 3 f* + b 4 f + ...)], where f is the frequency, k is 
wavenumber where k = 2*f/c and c is the velocity of light, r is an element-to-element distance, j is an imaginary 
unit, and ao, a 1f ... and b 0 , b v ... are the approximation coefficients. 

22 An apparatus according to claim 1 8, wherein the approximate equations for the mutual reactions By are B % = exp 
* {-ikr)[c 0 + Cl f2 + c 2 f 4 + c 3 f6 + c 4 f 8 + . .. + j{d 0 f + d,!* + d 2 f5 + d 3 f? + d 4 f* +...)]. where f is the frequency, k is wavenumber 

where k = 2«f/c and c is the velocity of light, r is an element-to-element distance, j is an imaginary unit, and Cq, 
c 1t ...and d 0 , d-,, ..'I are the approximation coefficients. 

23. An apparatus according to claim 15, wherein the number of sample frequencies is equal to the number of terms 
in a real part and in an imaginary part of each of the approximate equations. 

24. An apparatus according to claim 23, further comprising means for determining values of the plurality of sample 
frequencies based on a specified frequency range and on the number of sample frequencies. 

25. A program storage device readable by a machine, tangibly embodying a program of instructions executable by 
' the machine to perform method steps for calculating an electromagnetic field intensity, the method steps compris- 
ing: 

a) calculating approximation coefficients for each of approximate equations that respectively yield approximate 
values of electromagnetic characteristic values, between a plurality of elements constituting an electric circuit 
device and of the elements themselves, at an arbitrary frequency; and 

b) calculating the approximate values of the electromagnetic characteristic values between the elements and 
of the elements themselves at a designated frequency in accordance with the approximate equations having 

• the approximation coefficients calculated in step a), thereby making it possible to evaluate the electromagnetic 
field intensity for the designated frequency. 

26. A program storage device according to claim 25, wherein the method steps.further comprise the steps of: 

• c) calculating an electric current flowing on each element by solving simultaneous equations obtained by sub- 
stituting the approximate values of the electromagnetic characteristic values calculated in step b) into equations 
respectively describing boundary conditions for the elements; and 

d) calculating the electromagnetic field intensity for the designated frequency from the electric currents re- 
spectively flowing on the elements calculated in step c). 

27. A program storage device according to claim 25, wherein the step a) includes the substeps of: 

i) calculating the electromagnetic characteristic values between the elements and of the elements themselves 
from geometric data of the elements at each of a plurality of sample frequencies; and 

ii) calculating the approximation coefficients by solving simultaneous equations obtained by substituting the 
electromagnetic characteristic values calculated at the plurality of sample frequencies into each of the approx- 
imate equations. 

28. A program storage device according to claim 27, wherein the substep a) ii) includes the substeps of: 

calculating matrix elements in a lower triangular matrix and an upper triangular* matrix from the sample fre- 
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quencies, said lower triangular matrix and upper triangular matrix being derived from a matrix containing pow- 
ers of the sample frequencies as matrix elements by being decomposed as a product of a lower triangular 
matrix, with diagonal elements of 1 , and an upper triangular matrix; and 

calculating the approximation coefficient from the electromagnetic characteristic values and the matrix e IB- 
S' ments in the lower triangular matrix and the upper triangular matrix. 

29. A program storage device according to claim 25, wherein the electromagnetic characteristic values include mutual 
impedances between the elements constituting the electric circuit device and of the elements themselves. 

10 30. A program storage device according to claim 29, wherein 

the elements constituting the electric circuit device include conductor elements and dielectric elements, and 
the electromagnetic characteristic values further include mutual admittances between the dielectric elements 
and of the dielectric elements themselves, mutual reactions between the conductor and dielectric elements, 
is and mutual reactions between the dielectric elements and of the dielectric elements themselves. . 

31. A program storage device according to claim 30, wherein the method steps further comprise the steps of: 

c) calculating an electric current, an equivalent electric current, and an equivalent magnetic current flowing 
20 on each element, by solving simultaneous equations obtained by substituting the approximate values of the 

electromagnetic characteristic values calculated in step b) into equations respectively descrtoing boundary 
conditions for the elements; and 

d) calculating the electromagnetic field intensity for the designated frequency from the electric currents, the 
equivalent electric currents, and the equivalent magnetic currents respectively flowing on the elements, cal- 

25 culated in step c). 

32. A program storage device according to claim 29, wherein the approximate equations for the mutual impedances 
Zjj are Z t] = exp(-jkr)[a 0 + a,f 2 + a 2 f 4 + a 3 f6 + a 4 f* + . .. + j(b(/f + b,f + b 2 f 3 + b 3 f 5 + b 4 f 7 + ...)], where f is the frequency, 
k is wavenumber where k - 2^f/c and c is the velocity of light, r is an element-to-element distance, j is an imaginary 

30 unit, and ao, a v .... and b 0 , b v ... are the approximation coefficients. 

33. A program storage. device according to claim 30, wherein the approximate equations for the mutual admittances 
Yy are = exp(-jkr)[ao + a,f2 + aaf 4 + a 3 f 6 + a^ + ... + Kb^f + b^ + b 2 f 3 + b 3 f 5 + b 4 f 7 + ...)], where f is the frequency, 
k is wavenumber where k = 27tf/c and c is the velocity of light, r is an element-to-element distance, j is an imaginary 

35 unit, and a^ a v ... and b 0 , b v ... are the approximation coefficients. 

34. A program storage device according to claim 30, wherein the approximate equations for the mutual reactions By 
are B i} = exp(-jkr)[c 0 + c^ 2 + c 2 f 4 + c 3 f 6 + c 4 1* + ... + j(d 0 f + d^ 3 + d 2 f 5 + d 3 t 7 + d 4 f 9 + ...)]. where f is the frequency, 
k is wavenumber where k = 2itf/c and c is the velocity of light, r is an element-to-element distance, j is an imaginary 

40 unit, and c 0 , c v ... and d 0 , d v ... are the approximation coefficients. 

35. A program storage device according to claim 27, wherein the number of sample frequencies is equal to the number 
of terms in a real part and in an imaginary part of each of the approximate equations. 

45 36. A program storage device according to claim 35, wherein the method steps further comprise the step of determining 
values of the plurality of sample frequencies based on a specified frequency range and on the number of sample 
frequencies. 
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